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Resumen 


El objetivo del presente estudio fue analizar la influencia de diferentes 
indicadores fisicoquímicos y microbiológicos de estructuras de la 
comunidad bacteriana en un ecosistema de agua dulce ubicado en la 
Reserva de la Biosfera Sierra del Rosario. La composición de las 
comunidades de bacterioplancton se determinó por electroforesis en 
gradiente desnaturalizante (DGGE) del 165 rDNA durante tres campañas 
en tres estaciones de muestreo en el río San Juan. Se realizaron análisis 
estadísticos multivariados de los patrones de la electroforesis en relación 
con parámetros químicos, biológicos y físicos. Los análisis mostraron la 
variación estacional de la comunidad bacteriana en las estaciones de 
muestreo menos contaminadas. Además, los sólidos disueltos totales 
(TDS), la concentración de E. coli y el oxígeno disuelto influyeron 
significativamente en los perfiles de DGGE en el río San Juan, lo que 
indica que al menos dos de las tres estaciones de muestreo analizadas 
tuvieron una buena calidad del agua. Con la combinación de variables 
ambientales y la DGGE es posible analizar la influencia de la 
contaminación química y microbiana en la comunidad bacteriana del río 
San Juan. En los sitios menos contaminados, las bajas concentraciones 
de TDS y los indicadores bacterianos de contaminación se relacionaron 
con la diversidad bacteriana. La presente investigación constituye el 
primer análisisindependiente de cultivo de las comunidades bacterianas 
de ecosistemas dulceacuícolas en una reserva de la biosfera en las islas 
del Caribe. 
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Palabras clave: electroforesis en gel de gradiente desnaturalizante 
(DGGE), ecosistema de agua dulce, país tropical, análisis multivariado, 


reserva de la biosfera. 


Abstract 


The objective of the present study was to analyze the influence of 
different physicochemical and microbiological indicators on bacterial 
assemblagesin a freshwater ecosystem located in a Biosphere Reserve. 
The bacterioplankton communities composition was determined by 
denaturing gradient gel electrophoresis (DGGE) of 165 rDNA during 
three sampling campaigns in three sampling stations of San Juan River. 
Multivariate statistical analyses of gel patterns, in relation to chemical, 
biological and physical parameters were performed. The analyses 
showed a seasonal variation of bacterial community at the less polluted 
stations. The correspondence canonical analysis showed that total 
dissolved solids (TDS), E. coli concentration and dissolved oxygen were 
significantly influencing the DGGE profiles at San Juan River, indicating 
that at least two of the three sampling stations analyzed possessed good 
water quality. With the combination of environmental variables and the 
DGGE molecular data it could be appreciated the influence of chemical 
and microbial contamination on the bacterial community of San Juan 
River. At the less polluted stations low concentrations of TDS and 
bacterial indicators of fecal contamination were related with bacterial 
diversity. The present investigation constitutes the first culture- 
independent analysis of freshwater bacterial communities in a Biosphere 
Reserve of the Caribbean islands. 
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Introducción 


Las comunidades bacterianas son muy sensibles a los cambios 
ambientales, así como a los cambios persistentes en los ecosistemas 
debido a la contaminación de la actividad humana o al cambio climático 
global (Lear, Anderson, Smith, Boxen, € Lewis, 2008; Washington, Lear, 
Neale, € Lewis, 2013). Investigaciones recientes sobre comunidades 
bacterianas sugieren que pueden ofrecer una herramienta adicional para 
evaluar la calidad ecológica de los cursos de agua (Lear et a/., 2012; 
Washington et a/., 2013; Newton €: McLellan, 2015; Wang, Chen, Yuan, 
Li, 8 Li, 2016). Esto se debe a las ventajas que ofrecen las comunidades 
bacterianas como indicadores, ya que pueden experimentar cambios en 
su composición en respuesta a perturbaciones ambientales, siendo 
detectados los efectos de estos cambios antes de que se transmitan a 


niveles superiores de las cadenas tróficas en los ecosistemas acuáticos 
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(Lewis et al., 2010). Además, pueden muestrearse en sitios donde no 
están presentes otros bioindicadores, como invertebrados bentónicos y 
peces (Lear et al., 2012). 


Los métodos independientes del cultivo contribuyen a reducir el 
sesgo en el análisis de la riqueza de las comunidades bacterianas en 
comparación con los métodos dependientes del cultivo debido a la 
presión selectiva que imponen estos últimos para crecer en un medio de 
cultivo sólido, favoreciendo el crecimiento de aquellas poblaciones 
dentro la comunidad bacteriana capaz de formar colonias. Es por esta 
razón que los métodos dependientes del cultivo subestiman la riqueza 
de especies en una muestra, que se estima mejor a través de métodos 
independientes del cultivo que pueden detectar bacterias viables, pero 
no cultivables (Burtscher, Zibuschka, Mach, Lindner, € Farnleitner, 
2009; Vaz-Moreira, Egas, Nunes, 8 Manaia, 2011). 


Los métodos moleculares como la electroforesis en gel con 
gradiente desnaturalizante (DGGE), el polimorfismo de longitud del 
fragmento de restricción terminal (T-RFLP) y el análisis de espaciador 
intergénico ribosómico automatizado (ARISA) permiten establecer 
perfiles comunitarios y afiliaciones filogenéticas de las especies 
detectadas en caso de que los fragmentos de ADN se puedan 
secuenciar. Particularmente, la DGGE permite obtener información 
taxonómica, ya que las bandas obtenidas se pueden cortar, reamplificar 
y secuenciar. Así, esta técnica combina las ventajas de la clonación, 
secuenciación y T-RFLP, y puede utilizarse para el análisis de 
comunidades bacterianas de diferentes ecosistemas (Duarte, Cassio, € 
Pascoal, 2012; Liu, Huang, Sun, Xu, € Xu, 2012). 
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El sector hidromineral "Las Terrazas" está ubicado en el municipio 
de La Palma, provincia de Artemisa, a orillas del ío San Juan, dentro del 
área de desarrollo del Complejo Las Terrazas, en la parte noreste de la 
Reserva de la Biosfera Sierra del Rosario, Cuba. Las aguas del sector 
son representativas de la zona de saturación de los macizos 
carbonatados, y emergen en el contacto con fallas paralelas y 
transversales a la Falla Pinar (Peña, Fagundo, Delgado, € Orbera, 
2001). Teniendo en cuenta la mineralización global, las aguas se 
consideran mineromedicinales, y según su uso se clasifican como aguas 
de uso balneológico y terapéutico (Peña et al., 2001). Si bien la 
evaluación de la calidad química y microbiológica del río San Juan 
(Arpajón, Larrea, Rojas, Heydrich, 8 Lugo, 2012; Larrea, Rojas, Lugo, € 
Heydrich, 2013; Romeu, Quintero, Larrea, Rojas, 8 Heydrich, 2015) se 
ha realizado antes, nunca se han implementado métodos independientes 
del cultivo para evaluar la comunidad bacteriana de dicho ecosistema. 
Este tipo de investigación podría ser beneficioso para conocer desde el 
punto de vista sinérgico la influencia de diferentes variables 
fisicoquímicas en la estructura de las comunidades bacterianas, las 
cuales pueden contribuir a la mejora de la gestión de tal ecosistema de 
agua dulce, de gran importancia para el turismo y el medio ambiente. 
Este estudio tiene como objetivo analizar la influencia de diferentes 
indicadores fisicoquímicos y microbiológicos sobre las comunidades 
bacterianas en un ecosistema de agua dulce ubicado en una reserva de 
la biosfera. 
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Materiales y métodos 


Muestreo y datos ambientales 


Para analizarla estructura de la comunidad de bacterioplancton en el río 
San Juan se realizaron muestreos en tres ocasiones (noviembre de 
2011; febrero de 2012, y abril de 2012) en tres estaciones de muestreo: 
presa El Palmar (S-1) (22% 50' 42.35" N y 829 56' 23.33" W); presa San 
Juan (S-2) (22% 50' 46.01" N y 82% 56' 27,09" W), y Baños del San 
Juan (S-3) (220 49' 24.02" N y 8290 55' 35.08" W). 


Para la toma de muestras se utilizaron frascos de muestreo de 
plástico estériles (1 1) (APHA, 2012). Todas las muestras fueron 
recolectadas por triplicado. Se tomaron muestras de agua a un metro de 
la orilla y 15 cm de profundidad. Los parámetros de campo que 
consisten en conductividad eléctrica, pH, temperatura, sólidos totales 
disueltos (TDS) y oxígeno disuelto (DO) se midieron en cada ubicación 
de muestra con el uso de un multímetro (Mettler Toledo). Todas las 
muestras de agua se transportaron en hielo al laboratorio y se 


analizaron en 12 horas. 


Las submuestras de agua utilizadas para el análisis de nutrientes 
(PO43--P, NHa4*-N, NO2"-N y NOs3"-N) y la demanda química de oxígeno 


(DQO) se filtraron con filtros GF/F secados al horno (Whatman 
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International Ltd., Maidstone, Reino Unido) según los protocolos 
estándar franceses AFNOR (2009). La determinación de nutrientes 
inorgánicos y DQO se hizo mediante métodos colorimétricos. En 
específico, la DQO se determinó mediante el método de digestión en 
reactor de dicromato, para lo cual se utilizó la cubeta de ensayo de DQO 
(Spectroquant 1 152 0001, rango 0-1500 mg.li DQO). Para la 
determinación de amoniaco se utilizó el método de Nessler y para el 
fosfato se utilizó el método del ácido ascórbico azul. Todos estos 
métodos se describen en la Asociación Estadounidense de Salud Pública 
(APHA) (APHA, 2012). Los datos de precipitación se obtuvieron de la 
estación Ecológica en la provincia de Artemisa, Cuba. 


Determinación de coliformes termotolerantes, 
Escherichia coli y enterococos 


Los coliformes termotolerantes (TtC) se cuantificaron mediante filtración 
por membrana (filtros de nitrato de celulosa estériles de 0.45 um de 
tamaño de poro y 47 mm de diámetro, Sartorius) en agar lactosa con 
Tergitol (concentración final de 0.095 %o peso/volumen) y trifenilo 
2,3,5-cloruro de tetrazolio (TTC) (concentración final de 0.024 %o p/vol) 
(AFNOR, 2009), Escherichia coli (EC) y enterococos (Ent) se 


cuantificaron por el mismo método en agar Coliform Chromocult y agar 
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Enterococci Chromocult, respectivamente de acuerdo con Ouattara, 
Garcia-Armisen, Anzil, Brion y Servais (2014). Las bacterias se 
cultivaron por 24 h a 37 *C para EC y Ent, y por 24 ha 44 *C para TtC. 
Las colonias naranjas que producían un halo amarillo debajo de la 
membrana después de la incubación se consideraron colonias de 
coliformes termotolerantes; las colonias de color azul oscuro se 
consideraron colonias de E. colí y las colonias rojas se consideraron 
enterococos. Los recuentos se expresaron como unidades formadoras de 
colonias (UFC) por 100 ml de muestra. 


Análisis de las comunidades bacterianas del río San 
Juan 


Extracción de ADN y amplificación de los fragmentos de 
ADNr 16S 


La biomasa microbiana se colectó a partir de la columna de agua y se 
concentró por filtración. Se filtraron alícuotas de cada muestra (1 |) por 
triplicado a través de filtros de nitrato de celulosa de 5 um y 
posteriormente a través de filtros de nitrato de celulosa de 0.22 um. A 


partir de filtros de 0.22 um, la extracción de ADN se realizó mediante el 
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método SDS-fenol-cloroformo (Oh et a/., 2011). Todas las muestras de 
ADN se purificaron con DNeasy Blood and Tissue Kit (Qiagen). 


Para aumentar la sensibilidad y facilitar el análisis de DGGE, se 
realizó una técnica de reacción en cadena de la polimerasa (PCR) 
anidada. Durante la primera amplificación, el tubo de PCR contenía 25 ul 
de la mezcla de reacción que consiste en una concentración de 1 X de 
OneTaq Quick-Load 2 X Master Mix mezclado con tampón estándar 
(BioLabs, Nueva Inglaterra), una concentración final de 0.2 uM de los 
cebadores universales 27F (5 -AGAGTTTGATCMTGGCTCAG-3“) y 1492R 
(5“-TACCTTGTTACGACTT-3 “) para amplificar el gen ADNr 165 (Zhou et 
al., 2009); 1 ul de ADN molde (“ 70 ng), y agua MIilliQ estéril. Después 
de un paso de desnaturalización inicial de 5 min a 94 “C, las muestras 
se amplificaron durante 30 ciclos de 30 sa 94 *C,30sa55*Cy90Osa 
72 *C, y un paso de extensión final de 5 min a 72 *“C. 


A partir de los productos de la primera reacción de amplificación 
se realizaron dos reacciones de amplificación independientes: una con 
los cebadores 341FGC clamp (5 *- 
GOCCCGCCGCGCGCGGCGGGCGGGGCGGG 


GGCACGGGGGGCCTACGGGAGG-3") y 534R (5“-ATGTACGAGG-3 “) que 
amplifican la región variable V3 del ADNr 16S del dominio Bacteria 
(Muyzer, De-Waal, 8 Uitterinden, 1993); y el otro con los cebadores 
517FGC  (5'"-GTGCCAGCAGCCGCGG-3) (Gich, Schubert,  Bruns, 
Hoffelner, 8 Overmann, 2005) y 1061R (5'-CACGACACGAGCTGACGAC- 
3") (Degnan € Ochman, 2012), que amplifican las regiones variables 
V4-V6 del ADNr 165 del dominio de Bacteria. Ambas reacciones se 


realizaron en un volumen final de 50 pul de la mezcla de reacción que 
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contiene una concentración de 1 X Quick-OneTaq de Load Master Mix 
con 2 X tampón estándar (BioLabs, Nueva Inglaterra), 0.2 uM de los 
cebadores, 1 ul del producto amplificado y agua MIilliQ estéril. Después 
de un paso de desnaturalización inicial de 5 min a 95 *C, las muestras 
se amplificaron durante 30 ciclos de 1 min a 95 *C, 1 min a 53 *Cy 2 
min a 72 *C, y un paso de extensión final de 7 min a 72 *C. Todas las 
reacciones de PCR se realizaron en un termociclador (ThermalCycler C 
1000 TM BIO-RAD). 


DGGE 


La DGGE se realizó mediante el sistema Ingeny phor U (Rafer, España). 
Los productos de PCR obtenidos de la reacción con los cebadores 
341FGC y 534R se cargaron en un gel de poliacrilamida al 10 % (p/v) en 
tampón TAE 1X (Tris 20 mM, acetato 10 mM, EDTA 0.5 mM pH 8); 
mientras que los productos obtenidos de la reacción con los cebadores 
517FGC y 1061R se cargaron en un gel de poliacrilamida al 7 %. Se 
prepararon geles de poliacrilamida con un gradiente desnaturalizante de 
30 a 60 % (donde 100 % desnaturalizante contenía urea 7 M y 40 % 
formamida). Toda la electroforesis se realizó a un voltaje de 200 V 
durante 4 h y 30 min a 60 *C. Los geles se tiñeron con SYBR Gold 
(dilución 1:10 000) durante 30 min. Los geles teñidos se fotografiaron 
en un gabinete transiluminador UV, con una cámara de video (filtro de 


posiciones de gabinete de luz Alpha Innotech Multimage). 
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Análisis estadísticos 


Basado en la presencia (1) o ausencia (0) de bandas individuales en 
cada Carril, se construyó una matriz binaria utilizando el programa 
PyElph 1.4 (Pavel 8i Vasile, 2012), y se hizo un análisis comparativo de 
los perfiles DGGE con el software PRIMER 6 (PRIMER-E, Ltd., Reino 
Unido) (Clarke 8 Goney, 2006). Se realizó un análisis de conglomerados 
utilizando el coeficiente de similitud de Sorensen y la técnica de 
vinculación promedio mediante el método de grupos de pares no 
ponderados con media aritmética (UPGMA, siglas en inglés). 


Para ¡investigar las relaciones entre la comunidad de 
bacterioplancton y las variables ambientales en el río San Juan, se usó 
una aproximación de gradiente directo, es decir, un análisis de 
correspondencia canónica (CCA) utilizando el software CANOCO, versión 
4.5 para Windows (Ter-Braak 8 Smilauer, 2002). Este método fue 
elegido después de un análisisinicial por análisis de correspondencia sin 
tendencia (DCA), que reveló que los datos exhibían una respuesta 
unimodal a las variables ambientales. Los datos se importaron utilizando 
el programa WCanolImp dentro del paquete CANOCO. Luego se utilizó el 
programa CANOCO para realizar CCA con escalado de especies en 


distancias entre muestras para que las muestras y las variables 
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ambientales formaran un biplot. Se agregaron variables explicativas 
hasta que la adición de otras variables no contribuyó significativamente 
(p < 0.05) a una mejora sustancial del poder explicativo del modelo. 
Para evaluar estadísticamente la importancia del primer eje canónico y 
de todos los ejes canónicos juntos, se usó una prueba de permutación 
de Monte Carlo con 999 permutaciones no restringidas bajo el modelo 


reducido. 


Finalmente, para representar biplots, se utilizó el programa 
CANODRAW dentro del paquete CANOCO. Las correlaciones por pares de 
rango de Spearman entre los diferentes conjuntos de variables 
ambientales probadas mencionadas anteriormente se realizaron antes 


de la CCA para eliminar las variables covariantes. 


Cálculo del índice de calidad de agua 


El índice de calidad del agua (WQI, por sus siglas en inglés) se calculó 
de acuerdo con Yisa y Jimoh (2010). Primero se calculó la escala de 
calificación de calidad para cada parámetro gi utilizando esta expresión: 


qí = (Ci/Si) x 100 
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donde (Ci) es la concentración de cada parámetro en cada muestra de 
agua y (Si) es el valor estándar permisible para cada indicador dado por 
WHO (WHO, 2003). Luego, el peso relativo (Wi) se calculó mediante un 
valor inversamente proporcional al estándar recomendado (Si) del 
parámetro correspondiente: 


Wi = 1/Si 
Finalmente, se calcula el WQI general: 
WQI general = 2qiWi¡/2 Wi 


Según Ramakrishnaiah, Sadashivaiah y Ranganna (2009) hay 
cinco categorías para la calidad del agua: excelente calidad del agua (< 
50); buena calidad del agua (50-100); mala calidad del agua (100-200); 
muy mala calidad del agua (200-300), e inadecuada (> 300). 


Resultados 


Calidad fisicoquímica y microbiológica del agua 
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La calidad fisicoquímica y microbiológica del agua del río San Juan en el 
periodo 2011-2012 se muestra en la Tabla 1. Basado en su 
conductividad eléctrica, concentraciones de fosfato y nitrato, las 
estaciones de muestreo del río San Juan pueden considerarse 
mesotróficas, según la clasificación de Sigee (2005). Todas las 
estaciones de muestreo tuvieron buena calidad microbiológica de agua 
de acuerdo con las normas cubanas (Norma Cubana 22, 1999). Sin 
embargo, la estación Baños del San Juan (S-3) tuvo altos valores de 
concentraciones de nitrato y fosfato en comparación con las otras dos 
estaciones. Según el WQI, el agua de las estaciones S-1 y S-2 se puede 
clasificar como agua de excelente calidad (Tabla 2). Sin embargo, el 
agua de la estación S-3 está clasificada como inadecuada. La 
temperatura correlacionó negativamente con los TDS (r = -0.71, p < 
0.05) y el amonio (r = -0.85, p < 0.01). El pH correlacionó 
positivamente con el nitrato (r = 0.72, p < 0.05) y negativamente con 
la DQO (r = -0.84, p < 0.01). La conductividad eléctrica mostró una 
correlación positiva con los TDS (r = 0.97, p < 0.01) y una correlación 
negativa con las precipitaciones (r = -0.77, p < 0.05). El oxígeno 
disuelto correlacionó negativamente con las precipitaciones (r= -0.82, p 
< 0.01), y la concentración de nitrito correlacionó positivamente con la 


concentración de fosfatos (r = 0.81, p < 0.01). 


Tabla 1. Datos ambientales registrados de las estaciones de muestreo 
del río San Juan. 
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S Los datos se obtuvieron de la Estación Ecológica Sierra del Rosario, Artemisa. Se 
tuvieron en cuenta las precipitaciones de 10 días antes de la fecha de muestreo. 


TtC: coliformes termotolerantes; EC: Escherichia coli; Ent: enterococos; TDS: sólidos 
totales disueltos; DQO: demanda química de oxígeno; DO: oxígeno disuelto; S-1: 
Estación Presa El Palmar; S-2: Estación Presa San Juan; S-3: Estación Baños del San 


Juan. 


Tabla 2. Valores de WQI calculados para el río San Juan. 


Datos ambien- | Valores WQI general 
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DQO 10.00 975.00/1008.00/850.00| 97.50|100.80| 85.00 
(mg.L"*) 
_ MEN 
Ñ | | AIN 
Ñ | HEN 


NO3.N 10.00 60.90 | 112.60|151.35 
(mg.L"*) 


0.10 
0.10 
2.00 
0 2.00 55.30 68.60 |110.60| 140.00 


TDS: sólidos totales disueltos; DQO: demanda química de oxígeno; DO: oxígeno 


disuelto; S-1: Estación Presa El Palmar; S-2: Estación Presa San Juan; S-3: Estación 
Baños del San Juan. 


Análisis de las comunidades bacterianas del río San 
Juan mediante DGGE 


La Figura la muestra el patrón de bandas obtenido por la DGGE de las 
muestras del río San Juan. El número de bandas por muestra varió 
entre 15 (S-3 noviembre de 2011) y 26 (S-2 abril de 2012); además, es 
posible observar diferencias en los perfiles de bandas de cada estación 
de muestreo en cada uno de los meses analizados, así como diferencias 


en la intensidad de las bandas. 
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Figura 1. Análisis DGGE de fragmentos de ADNr 165 amplificados por 
PCR con los cebadores 341FGC y 534R. (A) Imagen del gel obtenida de 
muestras del río San Juan recolectadas en el periodo noviembre 2011- 
abril 2012. Gradiente desnaturalizante utilizado: 30-60 %. (B) Análisis 
de conglomerados por método UPGMA basado en patrones de bandas de 
muestras del río San Juan. La similitud se expresó como porcentaje del 
índice de Sorensen. Estaciones de muestreo del río San Juan: S-3: 
Baños del San Juan; S-2: Presa San Juan; S-1: Presa El Palmar. 


Teniendo en cuenta los perfiles de bandas obtenidos por DGGE, se 
realizó un análisis de agrupamiento (Figura 1b), donde se observó la 
formación de tres grupos con 60 % de similitud. El primer grupo 
compuesto por muestras de la estación S-3 recolectadas en febrero y 


abril de 2012; el segundo grupo incluye muestras de la estación de 
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muestreo S-1 recolectadas en los meses de febrero y abril de 2012, y el 
tercero agrupa las muestras de la estación S-2 recolectadas en febrero y 
abril de 2012. Se separaron las muestras del mes de noviembre de cada 
una de las estaciones analizadas. El análisis de permutación de varianza 
(PERMANOVA) indicó diferencias significativas entre las comunidades 
bacterianas de las diferentes estaciones de muestreo (p = 0.001), así 
como entre los diferentes meses muestreados (p = 0.001). Este 
resultado se puede ver en el análisis de conglomerados donde cada 
grupo formado está representado por una estación en particular; 
además, existe una separación entre las muestras recolectadas en los 
meses de febrero y abril de 2012, y las muestras obtenidas en 
noviembre de 2011. 


Análisis de la relación entre la estructura de la 
comunidad bacteriana, y los indicadores fisicoquímicos 
y de contaminación fecal 


Los sólidos totales disueltos (p = 0.015) y la concentración de E. coli (p 
= 0.045) fueron los principales factores que influyeron en la variabilidad 
de la estructura de las comunidades bacterianas en el río San Juan 
(Figura 2, Tabla 3). La prueba de Monte Carlo para todos los ejes 


canónicos (p = 0.002) fue significativa, lo cual indica que las variables 
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seleccionadas pueden explicar la variabilidad de la estructura de la 
comunidad bacteriana en las estaciones del río San Juan. Los sólidos 
totales disueltos, la concentración de E. coli y el oxígeno disuelto 
explicaron el 21.2, 19.4 y 17.3 %, respectivamente, de la variabilidad 
del bacterioplancton en el río San Juan. De acuerdo con el análisis de 
correspondencia canónica, el porcentaje acumulado de varianza 
explicado por las especies-variables ambientales indicó que el primer y 
segundo eje canónico describen el 76.8 % de la varianza. Las 
correlaciones especie-ambiente de cada eje fueron 0.93 y 0.98 para el 
primer y segundo eje, respectivamente. El oxígeno disuelto y los TDS 
fueron los dos componentes principales del primer eje, y la 
concentración de E. coli fue el componente principal del segundo eje. 
Basado en el patrón en el biplot, dividimos las muestras en tres grupos: 
Grupo 1, la mayoría de las muestras de las estaciones S-1 y S-2 (sin 
incluir noviembre de 2011) se caracterizaron por los valores más bajos 
de TDS y concentración de E. coli; Grupo 2, muestras de la estación S-3 
(de noviembre de 2011 a abril de 2012) se caracterizaron por los 
valores más altos de concentración de TDS y E. coli, y Grupo 3, que 
incluyó muestras de las estaciones S-1 y S-2 (noviembre de 2011), 
caracterizadas por los valores más altos de oxígeno disuelto. 
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Figura 2. Biplot de análisis de correspondencia canónica (CCA) de la 
estructura de las comunidades bacterianas del río San Juan en el 
periodo noviembre 2011-abril 2012 utilizando indicadores fisicoquímicos 
y microbiológicos. Las flechas rojas indican la dirección en la que 
aumenta el valor de cada indicador de calidad del agua. La longitud de 
las flechas indica el grado de correlación con el eje mostrado. La 
posición de las estaciones de muestreo en las flechas se interpreta 
mediante la proyección de los símbolos en las flechas e indica el grado 
de influencia del indicador de calidad del agua representado por la flecha 
roja en la estructura de la comunidad bacteriana de la estación de 


muestreo. DO: oxígeno disuelto; TDS: sólidos totales disueltos; EC: 
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Escherichia coli; S-1: Estación Presa El Palmar; S-2: Estación Presa San 
Juan; S-3: Estación Baños del San Juan. 


Tabla 3. Resumen de los resultados de los análisis de correspondencia 
canónica (CCA) de los datos de la estructura de la comunidad de 
bacterioplancton cuando están limitados por los indicadores de calidad 
del agua en el río San Juan. 
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Las letras en negrita indican los factores fuertemente correlacionados con los dos 
primeros ejes de ordenación de CCA (r > 0.70 o r > -0.70). 


Discusión 


Calidad fisicoquímica y microbiológica del río San Juan 


En el contexto de la evaluación del impacto de los contaminantes en la 
diversidad de comunidades microbianas en los ecosistemas de agua 
dulce, la diferencia entre los efectos relativos de la presión selectiva 
debido a la contaminación, y los que resultan de los factores y procesos 
ambientales naturales es muy importante. Se pueden encontrar 
variaciones naturales en la estructura y función de las comunidades no 
impactadas cuando se comparan diferentes ecosistemas, o diferentes 


estaciones de muestreo dentro del mismo ecosistema, como refieren 
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Dorigo, Lefranc, Leboulanger, Montuelle y Humbert (2009). Por esta 
razón, en el presente estudio se evalúan estaciones contaminadas y no 
contaminadas en el río San Juan, para analizar la influencia de los 


parámetros ambientales en la estructura de la comunidad bacteriana. 


El río San Juan está ubicado en una zona rural y particularmente 
en la Reserva de la Biosfera Sierra del Rosario, donde un Programa de 
Desarrollo Sostenible activo contribuye a la protección del medio 
ambiente (Larrea et al., 2013). Las aguas de este entomo son muy 
limpias, con abundantes peces, algunas macrófitas y buen olor, aunque 
la mayoría de los indicadores fisicoquímicos y microbianos no superan 
las normas cubanas para baños en cuerpos de agua costeros y 
continentales (Norma Cubana 22, 1999). La estación S-3 presentó los 
valores más altos de las variables microbianas y químicas, en 
comparación con las otras dos estaciones. Este hecho podría ser posible 
por el aporte antrópico que genera una gran afluencia de público a estas 
piscinas naturales y el efecto de la escorrentía de una laguna de 
oxidación, que se refleja en las concentraciones de amonio e indicadores 
fecales. Considerando el WQI en el periodo evaluado, las estaciones de 
muestreo S-1 y S-2 presentaron una excelente calidad de agua en 
comparación con la estación S-3, lo que está en correspondencia con 
estudios realizados previamente por Larrea et al. (2013). Para 
simplificar, informar e interpretar la información compleja obtenida de 
un cuerpo de agua, diferentes autores han calculado índices de calidad 
del agua (Yisa € Jimoh, 2010; Rubio-Arias, Contreras-Caraveo, 


Quintana, Saucedo-Teran, 8 Pinales-Munguia, 2012; Rubio-Anias et al., 
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2016). Así, con un solo número se puede interpretar el nivel de 
contaminación que presenta un cuerpo de agua. 


Basado en los coeficientes de correlación obtenidos en el presente 
estudio, la temperatura se relacionó negativamente con los sólidos 
totales disueltos y las concentraciones de amonio, lo que podría deberse 
a que el aumento de temperatura puede elevar la tasa metabólica de los 
microorganismos, lo que lleva a la degradación de la materia orgánica, 
que se utiliza como fuente de carbono y energía (Dinsmore, Billett, 8 
Dyson, 2013). Esto provoca la reducción de materia orgánica en la 
columna de agua y por lo tanto del total de sólidos disueltos. Por otro 
lado, el aumento de temperatura puede favorecer la actividad de las 
bacterias nitrificantes, que pueden utilizar el amonio en el proceso de 
nitrificación, que consta de dos pasos: (1) la oxidación del amonio a 
nitrito por las arqueas y bacterias oxidantes del amonio, y (2) oxidación 
de nitrito a nitrato llevada a cabo por bacterias oxidantes de nitrito 
(Chen, Cao, Song, 4 Zhou, 2010). De esta forma, con el aumento de 
temperatura, disminuirían las concentraciones de amonio en las aguas 


en las condiciones del río San Juan en el periodo analizado. 


En el caso de la conductividad eléctrica y los TDS, ambos son una 
medida indirecta de la salinidad de una solución (Niekerk, Silberbauer, € 
Maluleke, 2014), lo que explica la correlación positiva observada entre 
ambos indicadores en el río. La conductividad eléctrica se utiliza a 
menudo como un sustituto de los sólidos totales disueltos y es una 
medida de la capacidad de una solución para conducir corriente eléctrica 
entre dos puntos (Dunlop, McGregor, € Horrigan, 2005; Niekerk et a/., 


2014). Por el contrario, los TDS proporcionan una indicación general del 
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nivel de sólidos disueltos de menos de dos micrómetros presentes en 
una corriente fluvial, que incluye todos los electrolitos disociados que 
aumentan la concentración de sales y otros compuestos, como la 
materia orgánica disuelta (Tanor, Ts'enoli, 8 George, 2014) Además, se 
observó una correlación negativa entre la conductividad eléctrica y la 
precipitación, lo que puede deberse a que la conductividad eléctrica 
tiende a cambiar en función de la recarga del ecosistema: durante los 
períodos lluviosos tiende a disminuirsu valor debido a la dilución de las 
sales; mientras que durante los periodos secos ocurre lo contrario 
(Rashed € Younis, 2012). 


La DQO se relacionó negativamente con el pH, lo que puede 
deberse a que a mayor cantidad de materia orgánica en el ecosistema, 
mayor cantidad de oxígeno disuelto se necesita para descomponerlo, 
proceso que contribuye a la liberación de dióxido de carbono, que 
disminuye el pH del agua. Varol, Gókot, Bekleyen y Sen (2012) 
argumentan que la descomposición de la materia orgánica conduce a la 


acidificación y por lo tanto a una disminución del pH. 


En el río San Juan se observó una correlación negativa entre las 
concentraciones de oxígeno disuelto y las precipitaciones. Jung et al. 
(2014) obtuvieron resultados similares en el río Han, en Corea del Sur. 
Estos autores observaron que las precipitaciones se relacionaban 
negativamente con el oxígeno y positivamente con la concentración de 
sólidos en suspensión. Las precipitaciones pueden contribuir a aumentar 
la turbidez en los sistemas acuáticos porque provocan afloramientos de 
sedimentos a partir del fondo del río y arrastran los suelos adyacentes; 


esto disminuye la penetración de la luz solar, lo que afecta a la actividad 
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fitoplanctónica y por tanto al proceso de fotosíntesis, con la consiguiente 
disminución del oxígeno disuelto. Así, un aumento de la turbidez del 
agua puede provocar una disminución en el proceso de fotosíntesis, y 
con ello una disminución del pH y las concentraciones de oxígeno 


disuelto. 


Por otro lado, hubo una correlación positiva entre las 
concentraciones de nitrito y fosfato, que fueron mayores en la estación 
S-3, en comparación con las estaciones S-1 y S-2. Una causa que podría 
contribuir a estas concentraciones de nutrientes son los escurrimientos 
de la laguna de oxidación aeróbica que se ubica entre las estaciones S-2 
y S-3, que presenta roturas y recibe volúmenes de aguas residuales 
domésticas superiores a su capacidad de retención, como se observó 
durante el muestreo. En este mismo periodo se observó una correlación 
positiva entre las altas concentraciones de nitratos y el pH, resultados 
similares a los obtenidos por Zhang, Wang y Wu (2014) en el lago Erhai 
en China. Estos autores sugieren que a pH 8 se favorece el crecimiento 
y metabolismo de las bacterias nitrificantes. Además, en condiciones 
aeróbicas, la liberación de amoniaco de los sedimentos disminuye y éste 
se convierte en nitrato por nitrificación. Tal fenómeno contribuye a la 


reducción de amonio y al aumento de nitratos en las aguas. 


Relación entre la estructura de la comunidad 
bacteriana, y los indicadores fisicoquímicos y de 
contaminación fecal en el río San Juan 
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La combinación de DGGE con métodos multivariados permitió apreciar la 
variabilidad espacio-temporal en la estructura de las comunidades 
bacterianas durante el periodo de estudio, lo que se corresponde con 
resultados obtenidos por diferentes autores en distintos ecosistemas de 
agua dulce utilizando esta técnica molecular (Berdjeb, Ghiglione, € 
Jacquet, 2011a; Figueiredo et a/., 2012; Ibekwe, Leddy, Bold, € Graves, 
2012; Tang etal., 2014). 


Basado en la estructura de las comunidades bacterianas, se 
observó que el mes de noviembre de 2011 se separa de los meses de 
febrero y abril de 2012 (Figura 1b). Además, las estaciones menos 
contaminadas (S-1 y S-2) se separan de la estación con mayor 
contaminación (S-3), lo cual apoya la idea de que las comunidades 
bacterianas pueden utilizarse como indicadores del estado ecológico de 
los ecosistemas acuáticos (Lear et a/., 2012; Washington et a/., 2013). 
Se argumenta que las comunidades bacterianas son muy sensiblesa los 
cambios ambientales, así como a los cambios persistentes en los 
ecosistemas debido a la contaminación producto de la actividad humana 
o al cambio climático global (Lear et a/., 2008). Las modificaciones 
temporales y espaciales detectadas pueden estar relacionadas con 
factores abióticos, como temperatura, oxígeno disuelto, pH, 
conductividad eléctrica, transparencia del agua, concentración de 
materia orgánica y nutrientes, entre otros; así como con factores 
bióticos, como contenido de clorofila a, interacciones con fitoplancton y 


zooplancton, depredación, parasitismo y competencia, que regulan los 
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cambios comunitarios en los ecosistemas acuáticos (Lindstróm, Kamst- 
Van-Agterveld, £ Zwart, 2005; Berdjeb et al., 2011a; Berdjeb, Pollet, 
Domaizon, €: Jacquet, 2011b; Llirós et a/., 2014). 


Los factores que más influyeron sobre la estructura de las 
comunidades bacterianas del río San Juan en el periodo noviembre 
2011-abril 2012 fueron los TDS y las concentraciones de E. coli; esta 
última utilizada como indicador de contaminación fecal. Ambas variables 
presentaron valores altos en la estación S-3, en comparación con las 
estaciones S-1 y S-2, lo que influye en los cambios en el 
bacterioplancton. La presencia de altas concentraciones de TDS denotan 
el aporte de los suelos a través de los procesos naturales de arrastre y 
drenaje, pero también la resuspensión del sedimento, que ocurre con 
mayor frecuencia durante los periodos lluviosos (Crim, Schoonover, 8 
Lockaby, 2012), como sucedió durante los meses de febrero y abril de 
2012, donde se reportaron 86.0 mm y 124.9 mm de precipitación, 
respectivamente. El aporte del suelo puede contribuir no sólo al ingreso 
de bacterias al agua del río sino también al ingreso de compuestos 
orgánicos e inorgánicos, todo lo cual promueve cambios en la estructura 
de las comunidades bacterianas, como se observó en el río Brazos, 
Texas, EUA (Becker, Rodibaugh, Hahn, € Nowlin, 2016). Además, la 
resuspensión del sedimento en la columna de agua también contribuye a 
los cambios en la estructura de las comunidades bacterianas, pasando 
microorganismos sedimentarios y nutrientes a la columna de agua. Las 
altas concentraciones de TDS pueden contribuir a la disminución de la 
luz en la columna de agua y, por lo tanto, a la baja de la fotosíntesis. 


Esto provoca una merma en la producción de oxígeno disuelto, lo que 
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contribuye a cambios en la estructura de la comunidad bacteriana 
(Breton-Deval, Sanchez-Flores, Juárez, € Vera-Estrella, 2019). Por otro 
lado, altas concentraciones de TDS pueden provocar un aumento de la 
salinidad, lo que trae consigo cambios en las comunidades bióticas, una 
disminución de la diversidad y la exclusión de especies menos tolerantes 
a estas condiciones (Weber-Scannell 8 Duffy, 2007). En paisajes 
rurales, como la Sierra del Rosario, donde se ubica el río San Juan, las 
fuentes de contaminación fecal incluyen ganado, escorrentías agrícolas y 
filtraciones de estanques de oxidación; estos últimos utilizados en el 
tratamiento de aguas residuales (Crim et a/., 2012; Sun et al., 2016). 
Goto y Yan (2011) argumentan que el suelo es un importante reservorio 
de bacterias indicadoras de contaminación fecal, y que el intercambio 
entre el suelo y el agua durante las precipitaciones abundantes 
constituye un medio de transporte de bacterias fecales. Estas bacterias 
pueden llegaral suelo a partir de las deposiciones de la vida silvestre (p. 
ej., aves, ciervos, cerdos salvajes, etc.), así como de animales 
domésticos. Las fuentes de contaminación mencionadas anteriormente 
tienen un impacto en la calidad de las aguas de la estación S-3 y, por 
ende, en los cambios que ocurren en las comunidades bacterianas de 


esa estación de muestreo. 


Esta investigación es de gran importancia porque, por primera vez 
en Cuba, se realiza un análisis de la estructura de la comunidad 
bacteriana de las aguas de una reserva de la biosfera, como la Sierra del 
Rosario, y se relaciona con la calidad del agua, lo que contribuye a la 
conservación de dicho recurso y su entorno. Como parte de los 


esfuerzos ambientales requeridos para la preservación del agua como 
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recurso natural, es necesario contar con datos cuantitativos relacionados 
con la contaminación del agua (CITMA, 2009), por lo que el presente 
estudio constituye un aporte para la mejora del manejo de los 
ecosistemas acuáticos. El uso de comunidades bacterianas como 
indicadores de la calidad del agua es muy importante, pues permite la 
detección de cualquier alteración en el ecosistema en una etapa 
temprana antes de que se transmita a niveles superiores de la cadena 


trófica en los ecosistemas acuáticos (Lewis et a/., 2010). 


Conclusiones 


Con la combinación de variables ambientales y datos moleculares de la 
DGGE se pudo apreciar la influencia de la contaminación química y 
microbiana en la comunidad bacteriana del río San Juan. En las 
estaciones menos contaminadas, las bajas concentraciones de TDS e 
indicadores bacterianos de contaminación fecal se relacionaron con la 
diversidad bacteriana. Tales resultados contribuyen al diagnóstico de los 
ecosistemas acuáticos analizados y constituyen el primer paso para los 
planes de gestión del río. El análisis en este documento ha 
proporcionado información útil sobre el papel de ciertas variables 
ambientales que influyen sobre la estructura comunitaria en los ríos. 


Resultaíña de gran interés estudios adicionales que incluyan una 


Tecnología y ciencias del agua, ISSN 2007-2422, 12(3), 82-123. DOI: 10.24850/j-tyca-2021-03-03 


ES 2021, Instituto Mexicano de Tecnología del Agua 
Tecnología y Open Access bajo la licencia CC BY-NC-SA 4.0 (https: //creativecommons.org 


CienciaszAgua  ncsa/4.0/) 


identificación profunda de los taxones que componen la comunidad 
bacteriana en el ecosistema de agua dulce analizado. 
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